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ABSTRACT 

Direct imaging and characterization of extrasolar Earth-like planets is strongly im- 
pacted by the orbital inclination of the planet to be studied, as a combination of 
pure geometrical effects and the impact of exozodiacal dust. Here, we perform simu- 
lations to quantify the impact of a priori knowledge of inclination for the efficiency of 
a typical coronagraphic or occulter-based mission. The relative impact and comple- 
mentarity with prior knowledge of exozodiacal brightness down to achievable levels is 
examined and discussed. It is found that inclination has an even greater impact than 
the exozodiacal brightness, though the two have excellent complementarity. We also 
discuss different methods for inclination determination, and their respective applica- 
bility to the context of precursor science to an imaging mission. It is found that if 
technologically achievable, a combined effort to determine inclinations and exozodia- 
cal brightnesses with ground-based facilities would substantially increase the efficiency 
of a space-based dedicated mission to image and characterize Earth-like planets. 
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1 INTRODUCTION 

, The planar geometries of planetary systems are of funda- 
' mental importance for the detection and characterization of 
I planets. Radial velocity is biased towards systems that are 
■ close to edge-on, and the typically unknown inclination leads 
to an uncertainty in mass. Transiting systems can only be 
observed under very tight inclination boundaries. In direct 
imaging, a planet on a highly inclined orbit will spend at 
least part of its time inside of the inner working angle (IWA) 
of the working instrument. Gravitational microlensing detec- 
tions favour projected separations close to the Einstein ring 
radius, hence projection effects from inclination will play a 
role also in this case. By contrast, the detection efficiency 
of astrometry can be said to be relatively independent of 
inclination (e.g. Unwin et al. 2008). 

It follows from this importance of orbital inclination 
for detection with a wide range of methods, that an a pri- 
ori determination of this quantity could potentially lead to a 
significant boost in detection efficiencies of (otherwise) blind 
surveys. An indirect way of doing so would be to determine 
the rotational plane of the star. Since a star and its disk form 
primordially from the gravitational collapse of a rotating 
cloud with conservation of angular momentum, the plane of 
the disk should in general be similar to the rotational plane 
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of the star, and so the difference in the final inclination of 
planets formed in the disk and the stellar rotation should be 
small. Indeed, we know that this is the case in our own Solar 
system (e.g. Beck & Giles 2005). We also know from systems 
like e Eri (Saar & Osten 1997; Greaves et al. 2005) and Fo- 
malhaut (Le Bouquin et al. 2009) that the rotational plane 
of the star is consistent with the disk plane of the circum- 
stellar material. A projection of the mutual inclination can 
also be measured directly in transiting systems through the 
Rossiter-McLaughlin effect. The general population of stud- 
ied hot Jupiters so far is consistent with coplanarity (e.g. 
Fabrycky & Winn 2009). However, some systems do exist 
in which the planet is markedly non-coplanar (e.g. Narita 
et al. 2009, Winn et al. 2009). Probably, a violent dynam- 
ical history or Kozai influence from a stellar companion is 
responsible in these cases. 

A study by Beatty & Seager (2010) investigated the 
impact of a priori known inclinations for the study of tran- 
sits among nearby stars. Partly motivated by this study, we 
here investigate the impact of known inclinations for the pur- 
pose of direct imaging of Earth-like planets. Direct imaging 
strongly benefits from low inclinations (as opposed to tran- 
sits, where high inclinations are necessary). The reasons for 
this are twofold: Consider an edge-on system {i — 90°) and 
a face-on system {i — 0°), at the same distance and with 
an equal planet at equal semi-major axis on a circular orbit. 
The semi-major axis is such that, at maximum projected 
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separation, the planet is outside of the observer's IWA. In 
the face-on case, the planet is then permanently visible, but 
in the edge-on case, for some fraction of the time the planet 
is inside of the IWA and will remain undetected if observed 
at an unfortunate epoch. Furthermore, however, the differ- 
ence in detectability is emphasized to an even greater extent 
if the detection limit is set by exozodiacal dust. In the edge- 
on case, a larger portion of the exozodiacal emission is pro- 
jected into the line-of-sight such that at maximal projected 
separation, the detection limit is significantly worse than for 
the face-on case. The edge-on planet spends the majority of 
its time well within the maximal projected separation, where 
the exozodiacal emission is even yet worse. 

In this paper, we quantify these effects for typical obser- 
vational parameters of dedicated Earth-like exoplanet imag- 
ing missions in the visual range. The conclusions apply 
equally to occulter-based (e.g. Cash 2006; Janson 2007) and 
coronagraph-based (Guyon et al. 2006) designs, although 
some practical differences exist which we will point out as 
they occur. The inclination is also of importance for infrared 
interferometric designs (e.g. Cockell et al. 2009), but these 
will not be considered in this paper. In the following, we will 
describe the input parameters to the simulations and how 
the simulations were performed, followed by the results and 
discussion. Some viable and non-viable methods for deter- 
mining inclinations are discussed before the conclusion. 



2 EFFICIENCY SIMULATION 

We have performed simulations to test the impact of know- 
ing the inclination distributions a priori, relative to lacking 
this knowledge. The idea is to simulate the relative time it 
would take to detect an Earth-like planet against an exo- 
zodiacal background at different points of time, for various 
inclinations, exozodi brightnesses, and target distances, as 
well as the fractional time that the planet is within and out- 
side of the IWA, respectively. These quantities are combined 
to quantify a relative observation efficiency between differ- 
ent targets under different circumstances. Different target 
selections are then made, which are blind to the inclinations 
or to various degrees of exozodiacal emission, and are com- 
pared to each other and to the perfect-knowledge case to 
determine the impact of various degrees of knowledge on 
the total efficiency of the survey. 

2.1 Input parameters and distributions 

The only constraints put on the observing facility are an 
IWA of 60 mas outside of which a contrast of lO"^'' - suffi- 
cient for an Earth-like planet to dominate the residual stellar 
signal - is reached, and a FWHM of 20 mas, corresponding 
roughly to a 6-meter telescope at 600 nm. The most im- 
portant parameter of the target sample is the distance, and 
so for a distance distribution of our simulated target stars 
we choose all targets from HIPPARCOS (Perryman et al. 
1997) within 16.7 pc with OA < B - V < 1.0 to represent 
reasonably solar-like conditions. For simplicity, the stars are 
all assumed to have the properties of a IMsun star, although 
in reality, the colour selection spans 0.7-1.3Msu„. Here the 
outer distance of 16.7 pc is simply set by the fact that be- 
yond this distance, the angular separation corresponding to 



1 AU is smaller than the IWA. Every star is assumed to 
have exactly one Earth-like planet at a circular 1 AU or- 
bit with a star-planet contrast of 3.34 * 10"^*^ (e.g. Janson 
2007). The random distribution of inclinations is set as the 
arccosine of a uniform random distribution, to accurately 
describe the real distribution. The number of stars fulfilling 
the distance and colour selection criteria is 164, and we will 
assume that 60 of those would be selected for the mission. 
The optimal number of stars in an actual survey depends 
on several parameters, particularly including the fraction of 
systems that host habitable planets, tiq. If 7)0 is known a 
priori, it is possible to optimize the number of targets for a 
given mission concept (e.g. Beichman 2003). 

For the statistical brightness distribution of exozodia- 
cal dust disks, we use a modified single-population equation 
from Greaves & Wyatt (2010): 

Z{x) = ex-" (1) 

where x is the rank of the star such that among 100 
stars, the brightest disk has x = 0.01 and the faintest 
has X — 1.00, and c is a normalization constant such that 
Z{0.9) = 1 zodi. The latter normalization means that at 
the rank of the Sun, which is assumed to have rank 0.9, the 
brightness is 1 zodi. 

The original exponent in Greaves & Wyatt (2010) was 
Q — 1.08. However, this relation was derived from infrared 
excess at 70 /im, which probes cold dust in the outer disk, 
and hence a different region of the disk from what we are 
concerned with here. A better indicator for warm dust in 
the inner disk is excess at 24 fj,m, and hence this wavelength 
range is preferable for our purpose of estimating the bright- 
ness distribution of exozodi. Excess at 24 ^m is often harder 
to detect than at 70 /im because of the more demanding con- 
trast to the stellar photosphere, but some detections were 
made by Beichman et al. (2006) and Trilling et al. (2008), 
the samples that also form the basis for the 70 /im excesses 
used by Greaves & Wyatt (2010). In total, 13 systems ex- 
ist which have 24 /im excess with a confidence of >3cr in 
these samples. Recently, Koerner et al. (2010) presented 29 
systems with 24 /xm excess, but these are primarily around 
K-type stars, whereas we are concerned with Sun-like stars. 
The Koerner et al. (2010) sample lacks systems of the or- 
der of the highest excesses in Beichman et al. (2006) and 
Trilling et al. (2008), which should be those that are easi- 
est to detect. Hence, we consider that there is a risk that 
the low-mass sample will bias the result and do not use this 
sample for the main simulation. We plot the 13 systems in 
Fig-U along with the a — 1.08 power law. The fit is reason- 
able, but cannot match the two highest excess cases. The 
full ensemble of stars can be fit better by using a somewhat 
shallower slope, the best fit is provided by a = 0.80, also 
plotted in the figure. Hence, we use this exponent in the 
simulations, although it is based on fewer targets that the 
relation fit in Greaves & Wyatt (2010), since 24 fim excess is 
more reasonable for the exozodiacal dust we are concerned 
with. 

It should be emphasized that even the original distri- 
bution is an extrapolation based on observations in a very 
limited parameter range. Hence, the distribution is very un- 
certain, but represents a "best guess" case on the basis of 
the existing observations. Note that alternative simulations 
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Figure 1. Infrared excess emission at 24 /im for known targets 
(stars) in the literature. The solid line corresponds to the 
power law derived from 70 /im excesses. A best-fit to the full 
ensemble (dashed line) is provided by x~'^'^^ (see the text for 
discussion) . 



based both on the 70 /tm excesses and the Koerner et al. 
(2010) 24 /xm excesses are presented in Sect. 12.41 

The face-on surface density distribution of each exozo- 
diacal disk is another uncertainty, but here it is assumed 
to be r~^'^ out to 3 AU. Studies in the visual wavelength 
range have implied that this is the case for the zodiacal disk 
at 0.3-1 AU, assuming a thin wedge-like structure (Leinert 
et al. 1981), and this holds approximately true also at ~1- 
3 AU (Manner et al. 1976). The reflected emission goes as 
r~^", so in total the brightness distribution goes as r~^"^ . 
However, a recent study at 24 /im has given a brightness 
distribution of about r~'^ (Nesvorny et al. 2010). Assum- 
ing that the dust temperatures are much higher than 120 K 
such that the radiation is in the Ray leigh- Jeans regime, the 
temperature of the dust emission will scale as r~°'*, which 
leads to a density distribution of r"^'^, which in turn leads 
to a visual brightness distribution of r"*'^, vastly different 
from what we have assumed above. Since we are primarily 
interested in visual rather than infrared radiation here, we 
adopt the former value for the main simulation. However, 
to emphasize the effect of this large discrepancy, we also in- 
clude a simulation based on the latter assumption in Sect. 
[2:41 

For normalization, the surface brightness of the zodiacal 
cloud at lAU (1 zodi) is 23 mag (arcsec)"'^ (e.g. Absil et 
al. 2009). The disk is assumed to be azimuthally flat. In 
reality, the dynamical influence of planets in the disk will 
introduce asymmetric features. A simulation based on the 
Earth case has been performed by Dermott et al. (1994). In 
those simulations, it was found that the excess flux from the 
strongest morphological feature (a density enhancement in 
an Earth-trailing orbit) is only 10%. Since our simulations 
span about a factor 1000 in disk brightnesses, this effect is 
negligible in comparison, and can be disregarded. The effect 
that it would have, if stronger, would be maximal for a disk 
seen edge-on, and minimal for face-on. 

The contribution of zodiacal dust from the Solar system 
itself within the telescope beam (/z) is assumed to be 25% of 
the brightness of a 1 zodi edge-on exozodiacal disk, following 
the values used in Agol (2007). 
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Figure 2. Flux fraction between an Earth-like planet and an exo- 
zodiacal disk as function of inclination, at distances of 5 pc (solid 
line), 10 pc (dashed line), and 15 pc (dotted line), at a projected 
separation of 1 AU. Even though the fraction is a strong function 
of distance, it is still the case that a very small inclination system 
at 15 pc gives a better sensitivity than a very high inclination 5 
pc system, even at the maximum projected separation (1 AU). 
The inclination is plotted in a cos i scale to represent uniform 
probabilities. The edge-on case is cosi = and the face-on case 
cos 1 = 1. The flattening close to edge-on is due to the 3 AU outer 
cut-off of the exozodiacal disk. 

2.2 Procedure 

For 1000 populations p and for 100 instances of time t per 
population, we generate random inclinations i and total ex- 
ozodiacal brightnesses z for the 164 targets according to the 
distributions described above. At each instance, the planet 
is located at a projected separation 

p{t) = yJ{cos(t){t)Y + (sin0(t) cosi)2 (2) 

where 4>{i) is the phase angle at some instance of the 
orbit. The projected separation is either outside of or inside 
of the IWA, setting /vis,t to one or zero, respectively. The 
average over instances /vis is the fractional time of an orbit 
during which the planet is visible. 

For each instance when the planet is visible, the line- 
of-sight brightness of the exozodiacal dust /ez,t must be cal- 
culated at the position of the planet. This is done by first 
calculating the surface brightness for a face-on disk as a 
function of physical position in two dimensions, in steps of 
0.002 AU out to 3 AU. The flux is then evaluated in a box 
corresponding to a 20 by 20 mas area. The box is centered 
on the X-axis, at a separation given by the projected sepa- 
ration of the planet at a given instance. A vertical stretch is 
also applied to the box to account for the tilt of the disk rel- 
ative to the observer. In other words, scaling as function of 
distance and projection as function of inclination are done 
by the corresponding inverse scaling and de-projection of 
the box. The flux is normalized by the randomly assigned z 
factor, and is finally normalized by . 

The planetary brightness /pi,t is calculated from the 
3.34* 10"^" contrast to an absolute magnitude of a Sun-like 
star of A4y = 4.7 mag, normalized to the distance of the 
system. Furthermore, we factor in an illumination function 
modeled by I{t) = 1 -)- sin(<^(i;)) sin(i) to account for the 
different phases of the planet during an orbit. Given both 
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the exozodiacal and planet brightness, we can calculate a 
signal-to-noise ratio as 



S_ 



pi,t 



-\//pl,i + /cz.t + /z + /ps 



(3) 



Here, /psf is the PSF brightness outside the IWA, which 
is 10"^'' of the stellar flux, as discussed above. From this we 
can calculate a relative time rt = (S/N)^^ required to reach 
some given sensitivity. With this parameter in hand, we can 
calculate a total inverse observing efhciency as ^ = r//vis 
for each target in each population. To summarize, the pa- 
rameter 5 corresponds to how long (relative to other targets) 
we would need to integrate to achieve an acceptable signif- 
icance for detection, with the inverse probability that the 
target will be observable in the first place (i.e., outside of 
the IWA) factored in. Note that we have chosen the median 
rather than the mean of r, this is because in nearby high- 
inclination systems, the rare situation can arise that a planet 
is formally visible (outside of the IWA), but the phase angle 
is such that hardly any flux reaches the observer, which leads 
to an absurdly long required integration time. When taking 
the mean, this rare instance affects all other instances to a 
dubious extent. Although it is formally correct under the 
assumption that the observer keeps integrating until find- 
ing an Earth-equivalent planet, or otherwise integrates long 
enough to have adequate sensitivity to detect the faintest 
Earth-equivalent object in the field, what would more re- 
alistically happen in these cases is that the observer would 
cut losses by ending the observation prematurely. Hence, 
the median is a more robust estimator for our purposes, al- 
though it may somewhat under-estimate the positive impact 
of prior knowledge. 

Using 5, we can select the 60 best targets for a survey. 
If we had perfect knowledge of the inclination and the exo- 
zodiacal brightness, we could make the perfect selection to 
minimize E, the sum of ^ for the selected 60 targets in ev- 
ery randomly generated population. We denote the average 
E resulting from a target selection with perfect knowledge 
over all populations as Sidoai- 

The relative efficiency resulting from different degrees 
of actual knowledge can now be evaluated by systematically 
losing information and calculating the E that results from 
the corresponding target selections. We assume that there 
are three degrees of knowledge for z. The worst knowledge 
is if we only know about all systems with z ^ 10^ zodi. This 
roughly corresponds to our present knowledge, and is de- 
noted as case 'z3'. The mid- level case of knowledge is when 
all systems with z J5 10^ zodi are known, denoted 'z2'. This 
represents a reasonable goal for a near-term dedicated effort 
to characterize exozodiacal disks. Case 'zl' with all z J5 10^ 
zodi systems known is the best case, and would require a 
highly substantial effort, but could be possible using the 
Large Binocular Telescope Interferometer (according to Ab- 
sil et al. 2008). For inclination, we assume that either no 
inclinations are known (case 'iO'), or that all inclinations 
are known (case 'if'). This is also a reasonable distinction, 
since except for a few special cases, the observational case for 
determining the inclination (discussed in Sect. [S]) should be 
about equal for nearby Sun-like stars, such that if the incli- 
nation can be determined in one case, it can be determined 
in all cases. 
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Figure 3. Dependence of ^ on distance d. The gray dots corre- 
spond to all targets in all populations of our simulation. Upper: 
The black symbols are the targets of one random population. The 
few systems that exist within ~5 pc are exceptional and will yield 
a low 5 almost regardless of z and i. Outside of this, the other 
parameters start to strongly influence ^. Lower: The dashed line 
shows the dependency on d for a hypothetical target with z = 1.6 
zodi and i = 60°. 



A pseudo-merit "if corresponding to the real merit E is 
calculated for all cases of limited knowledge. This is done by 
redoing the simulations and replacing every unknown quan- 
tity with the expected mean value of the statistical distribu- 
tion. For instance, in calculating \&z2,io, every i is replaced 
with 60°, and every z < 100 zodi is replaced with the value 
Z{0.5) — 1.6 zodi. A selection of 60 targets per popula- 
tion is done to minimize 5*22, io, and then the Sz2,io of the 
real distributions is calculated on the basis of that selection. 
In this way, the S of each case represents the real relative 
inverse efficiency based on a best-effort selection using the 
knowledge at hand. As a final step, we define the relative 
efficiency ecasc = Ez3,io/Ecasc. The normalization is done so 
that Ezs.io = 1-00 is the present-day efficiency. 



2.3 Results 

The values of each tcaao are listed in Table [T] A partic- 
ularly notable fact is that ez3,ii > ezi.io, i-e., knowing 
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Figure 4. Dependence of § on exozodiacal brightness z. The gray 
dots correspond to all targets in all populations of our simulation. 
Upper: The black symbols are the targets of one random popula- 
tion. This parameter is particularly important above z ~10 zodi, 
below this limit the g dependence is relatively weak. Lower: The 
dashed line shows the dependency on z for a hypothetical target 
with d = 10 pc and i = 60°. 



the inclinations has a larger impact on the mission effi- 
ciency than determining exozodiacal dust brightnesses to 
an achievable degree. The very best efficiency is achieved 
by doing both: ezi,ii is even close to eidoai- The fact that 
(ez3,ii — ez3,io) + (szi.io — Czs.io) is close to (ezi,ii - tzs.io) sig- 
nifies that there is considerable complementarity between 
the two factors, in the sense that one does not replace the 
other, but a maximal effort in both areas is what gives the 
optimal performance. Also notable is the fact that ezs.ii is 
23% higher than ez3,io. Since a program for determining the 
inclinations (if performed e.g. in the form of SONG, see Sect. 
13. 2|) is very likely to have significantly less than 20% of the 
cost of a dedicated mission searching for extrasolar Earth- 
like planets, it would clearly be a good investment to execute 
such a program in this context. 

The primary dependency on distance is due to the fact 
that within a fixed spatial sampling, the planetary flux de- 
creases with d^, whereas the disk brightness remains con- 
stant. For this reason, the few targets inside of 5 pc are ex- 
ceptionally good and will almost never decrease the average 
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Figure 5. Dependence of § on inclination i. The gray dots corre- 
spond to all targets in all populations of our simulation. Upper: 
The black symbols are the targets of one random population. 
The inclination is plotted in a cosi scale to represent uniform 
probabilities. The edge-on case is cosi = and the face-on case 
cosi = 1. There is a gradual improvement of ^ over the whole 
range from edge-on to face-on systems. Lower: The dashed line 
shows the dependency on i for a hypothetical target with d = 10 
pc and z = 1.6 zodi. 



Table 1. Relative efficiencies for various degrees of knowledge. 



z knowledge 


i unknown 


i known 


z > 1000 


ez3,i0 = 1-00 


ez3,il = 1-23 


2 > 100 


<:z2,i0 = 1-06 


ez2,il = 1-28 


2 > 10 


ezi,iO = 1-12 


^zi.il = 1-32 


Ideal case 


N/A 


^idcal = 1-33 



efficiency, regardless of the disk brightness and inclination. 
However, at 5 pc, the disk-to-planet contrast is already less 
favorable in the edge-on case than at 15 pc in the face-on 
case, even at the maximum projected separation of 1 AU. 
This is illustrated in Fig.[2l In the edge-on case, the contrast 
becomes yet worse as the projected separation decreases. 
Breakdowns of the dependency of the inverse efficiency ^ on 
d, z, and i are illustrated in Figs.OU andO respectively. As 
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mentioned, the distance is particularly important at small 
d. The exozodiacal brightness has the largest impact on ^ 
when z > 10. Below this, the dependency is much weaker, 
hence why ezi,ii is so close to eideai- It is worth noting that 
for the very most nearby stars, even a quite large z still has 
a quite small impact. The trail of data points at the bot- 
tom of the figure are due to different representations of the 
two most nearby stars (a Cen A/B). It can be seen that 
at the distance of 1.3 pc, the disk is so well-resolved that 
even at 2 ~ 100, the disk brightness is still dominated by 
the planet brightness, and does not contribute significantly 
to the noise. The inclination dependency extends over the 
whole range of i. 

In this simulation, we have chosen a set of realistic pa- 
rameters to provide an illustration of the effects of know- 
ing the inclination and exozodiacal brightness. The exact 
numbers can obviously be expected to vary somewhat with 
changed mission parameters or underlying physical assump- 
tions. Two important aspects - the residual PSF level and 
the exozodiacal brightness distribution - are discussed in 
the next section. Other mission parameters that may vary 
are the IWA or number of selected targets. The specific mis- 
sion parameters also affect how /vis should be factored in. 
For instance, in the event of a coronagraphic mission and 
with a priori known phase angles (e.g. from SIM, see Sect. 
I3.2|l . it is possible to schedule the observation such that it 
occurs when the planet is known to be outside of the IWA, 
and in this case, the present factoring of /vis is an over- 
correction. In the other extreme, if we have an occulter mis- 
sion where the target-to-target motion dominates the fuel 
budget such that the target list and visit schedule must be 
strictly pre-defined, then it becomes relatively more impor- 
tant that there is a high chance that the target is visible 
than what the required integration time is going to be, and 
then the present factoring of /vis is an under-correction. Fi- 
nally, specific target properties, such as binarity, exozodiacal 
dependence on mass etc. could also have some impact on the 
outcome. In any case, it is clear that the impacts of i and 
z are large, and that it is therefore highly useful to deter- 
mine these quantities a priori, to the greatest extent that is 
technologically feasible. 

2.4 Alternative results 

In the main simulation, we considered missions in which a 
contrast of 10"^" is reached outside of the IWA. However, 
since such a high contrast is technologically challenging, the 
idea is sometimes considered to loosen the constraints on 
the contrast of the residual stellar PSF, remove PSF sub- 
structure through differential imaging, and integrate against 
the residual random noise. Hence, we have also performed 
a calculation using a contrast of 10~® for comparison. The 
simulation was otherwise run with identical parameters as 
the main simulation. 

The results are shown in Table (2] where the relative 
efficiencies are denoted with index 'hp' to signify the high 
PSF case. Unsurprisingly, it can be seen that the impact of a 
priori knowledge decreases, though it remains significant, at 
least for the inclination. The decreased impact is due to the 
fact that the noise is to a much greater extent dominated by 
the residual stellar PSF, which depends on neither i nor z. 
However, the total efficiency of the mission also goes down 



Table 2. Alternative relative efficiencies in the event of high PSF 
residuals. 



z knowledge i unknown i known 



z > 1000 


•ihp.zS.iO = 


= 1.00 


fhp,z3,il = 


= 1.14 


z > 100 


<:hp,z2,i0 = 


= 1.02 


eiip,z2,il = 


= 1.16 


2 > 10 


^hp,zl,iO = 


= 1.04 


Elip,zl,il = 


= 1.17 


Ideal case 


N/A 




^lip, ideal - 


= 1.17 



Table 3. Alternative relative efficiencies in the event of a steep 
slope for the exozodiacal brightness distribution. 



2 knowledge 


i unknown 


i known 




2 > 1000 


est,z3,i0 = 1-00 


est,z3,il = 


1.33 


2 > 100 


fst,z2,i0 = 1-23 


Est,z2,il = 


1.59 


2 > 10 


est,zi,iO = 1-46 


^st.zl.il = 


1.77 


Ideal case 


N/A 


^st, ideal 


1.81 



dramatically, since what this scenario does is raise the noise 
fioor, and thus required integration time, for all targets. Due 
to the consequential rapid decrease of science-to-cost ratio 
with increasing residual PSF noise, it is dubious whether the 
high PSF case can ever be worthwhile for the purpose of a 
dedicated mission for finding Earth-like planets. 

The most uncertain physical conditions in the simula- 
tions are the brightness distributions of exozodiacal disks, 
both in terms of the z distribution and in terms of the 
brightness as function of r in individual disks. For the z 
distribution in the main simulations, we used a slope based 
on known 24 /im excesses, because it is a better representa- 
tive of warm dust in the inner disk than 70 /xm excesses. On 
the other hand, the number of Sun-like stars with known 70 
^m excesses is larger, since it is an easier quantity to mea- 
sure. The 70 fim excesses give a steeper slope of a = 1.08 
than the a = 0.80 slope derived at 24 ^m. For completeness, 
and to quantify the impact of this uncertainty, we run a set 
of simulations for the steeper slope case, denoted with index 
'st'. The resulting efficiencies are shown in Table [3] 

In this case, the impact of prior knowledge goes up, this 
is mainly because there are simply a larger number of stars 
with high values of z for the observer to identify and avoid. 
This also means that the relative importance of knowing z 
compared to knowing i goes up, for instance est,zi,io is higher 
than est,z3,ii. As mentioned above, we prefer the 24 /im slope 
for physical reasons, but on the other hand, it would be a 
disappointing experience if a large-scale dedicated mission 
for Earth-like planets was launched and it was found that the 
number of bright exozodiacal disks is larger than expected, 
so from a conservative viewpoint there are reasons to also 
consider the worst-case scenario. 

We also chose to use the Beichman et al. (2006) and 
Trilling et al. (2008) 24 jj,m excesses rather than those of 
Koerner et al. (2010), in order to avoid biasing the sample 
with the large number of K-type stars. On the other hand, 
the total Koerner et al. (2010) sample is larger, so in the 
event that it is representative, it is a better statistical in- 
dicator. Hence, we have also performed simulations with a 
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Table 4. Alternative relative efficiencies in the event of a shallow 
slope for the exozodiacal brightness distribution. 



z knowledge i unknown i known 



z > 1000 




iO 


= 1.00 


^sh.zS.il 


= 1.15 


2 > 100 


<=sh,z2,: 


iO 


= 1.00 


^sh,z2,il 


= 1.15 


2 > 10 


^sh.zl,: 


iO 


= 1.01 


Esh,zl,il 


= 1.15 


Ideal case 


N/A 






^sh, ideal 


= 1.16 



Table 5. Alternative relative efficiencies in the case where the 
morphology of the zodiacal disk is based on the thermal brightness 
distribution. 



z knowledge i unknown i known 



z > 1000 


fth,z3, 


iO 


= 1.00 


Eth,z3,il 


= 1.28 


z > 100 


<5th,z2, 


iO 


= 1.28 


Eth,z2,il 


= 1.59 


2 > 10 


^th.zl, 


iO 


= 1.58 


Eth.zl.il 


= 1.86 


Ideal case 


N/A 






^th. ideal 


= 1.92 



slope oi a — 0.36, which is the best fit to the Koerner et al. 
(2010) distribution. The results are shown in Table 21 

As expected, the total benefit of increased knowledge 
decreases in this case and reaches values similar to those of 
the high PSF case. The inclination remains an important 
parameter still, but the benefit of knowing z becomes prac- 
tically negligible. Hence, it is clear that the uncertainty of 
the exact z distribution leads to a broad possible range of 
mission impact of prior knowledge. From a risk minimization 
perspective, this in fact further underlines the importance 
of determining z and i a priori. 

Regarding the distribution of the disk brightness as 
function of r, we chose a slope of r~'^'^, on the basis of visual 
light observations. On the other hand, more recent mapping 
at 24 fim is rather more indicative of a slope of order r"'*'^, 
implying that our knowledge of this distribution is incom- 
plete. To take also this uncertainty into account, we perform 
yet another simulation with the r~*'^ slope. The numbers 
are provided in Table (5] 

Since the distribution is much steeper in this case, the 
benefit of prior knowledge increases significantly, providing 
the highest numbers among the simulations performed here. 
Given the fact that the brightness distribution of our own 
zodiacal disk is not only poorly known, but in addition may 
not necessarily be representative of the typical exozodiacal 
disk in the first place, it is necessary to consider that the 
typical distribution could be as steep as considered here, 
and thus the impact of prior i and z knowledge could be as 
large. 



3 INCLINATION DETERMINATIONS 

The fact that the inclinations of stellar systems are poorly 
known certainly indicates that such a determination is non- 
trivial. Beatty & Seager (2010) discuss some different meth- 
ods of inclination determination, here we will briefly discuss 
those for the context of a direct imaging survey as examined 



here, as well as discuss some additional methods. We sep- 
arate the discussion between small-number methods, which 
can be used to determine the inclination in special cases, 
and large-number methods, which could potentially deter- 
mine inclinations for the whole sample. The methods typ- 
ically assume that the stellar rotation and planetary orbit 
are coplanar. 



3.1 Small-number methods 

One method mentioned in Beatty & Seager (2010) is t;sini 
and Prot determinations to constrain sini. This is useful for 
particularly rapidly rotating stars, but probably not for typ- 
ical stars due to an ambiguity with macroturbulence. The 
analysis by e.g. Saar & Osten (1997) shows that even though 
the two effects are in principle distinguishable to some ex- 
tent, the resulting sini estimation is often of limited util- 
ity. A different method that could work in individual cases 
is determining the orbital plane of other circumstellar ma- 
terial, when such is present. We have already mentioned, 
e.g., the case of e Eri where a spatially resolved dust disk is 
present (Greaves et al. 2005). Outer giant planet compan- 
ions could be used for the same purpose. Such companions 
will be increasingly detectable over the coming years, with 
instruments like GPI or SPHERE, and in a longer perspec- 
tive, JWST and ground-based ELT:s. 

In principle, one might similarly expect that the orbital 
planes of multiple stellar systems could correspond to their 
rotational planes, and therefore the orbital planes of their 
respective planets. However, the formation and dynamical 
evolution of such systems leave uncertainties regarding the 
accuracy of inferring a rotational inclination from stellar or- 
bital inclinations. We have tested this accuracy on a set of 
known multiple systems. As an input sample, we use the 
Sixth Orbital Catalog of Visual Binary Star^j to select all 
stellar systems with three or more components. Out of these, 
we select only those systems where inclination and ascending 
node measurements exist with formal errors given, for two or 
more of the orbits in the system. One of the remaining sys- 
tems is Sgr A, which we remove as it is clearly not applicable 
to our study. What remains is 10 systems, of which three 
are quadruple (HIP 21402, Lane et al. 2007; HIP 28614, 
Muterspaugh et al. 2008; HIP 76563, Drummond et al. 1995) 
and the rest are triple systems (HIP 2552, Docobo et al. 
2008; HIP 9500, Hartkopf et al. 2000; HIP 79607, Raghavan 
et al. 2009; HIP 84949 and HIP 107354, Muterspaugh et al. 
2008; HIP 113726, Hartkopf et al. 2009; HIP 116726, Ole- 
vic & Cvetkovic 2005). The idea is to test whether at least 
the orbits within the system are coplanar with each other, 
which is a necessary pre-condition for a rotational plane to 
be inferred from any given orbit orientation. 

The differential orientation is calculated on the basis 
of spherical trigonometry from the differential inclination Ai 
and ascending node Afl: 

AC = arccos(cos(Ai) cos(Af2)) (4) 
For the quadruple systems (except HIP 76563 where 



^ http:/ /ad. usno.navy.mil/wds/orb6. html 
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Table 6. Orientations of triple and quadruple systems with mul- 
tiple known i and 







(Orbit 1 / Orbit 2) 


[deg] 


HIP 2552 AaAb/AB 


3.6±1.5 


HIP 9500 AB/AC 


50.5±2.5 


HIP 21402 AalAa2/AaAb 


125.9±2.0 


HIP 21402 AblAb2/AaAb 


106.3±5.8 


HIP 28614 AaAb/AB 


126.1±7.6 


HIP 28614 BaBb/AB 


93.5±3.8 


HIP 76563 AB/CD 


66.1±3.0 


HIP 79607 AB/BaBb 


175.0±3.2 


HIP 84949 AB/BaBb 


26.8±1.7 


HIP 107354 AB/BaBb 


47.8±2.5 


HIP 113726 AB/AaAb 


102.6±2.3 


HIP 116726 AB/BaBb 


28.7±9.5 



only two of the orbits are sufficiently constrained), two dif- 
ferential orientations are calculated, where the most pre- 
cisely determined orbit is related to each of the other two 
orbits. The results are shown in Table [6l It is immediately 
obvious that there is no strong preference towards copla- 
narity, although one system might be said to be conspicu- 
ously close to prograde coplanarity, and another one to ret- 
rograde coplanarity (but neither is perfectly coplanar within 
the errors). We test this more rigorously by generating 10^ 
orbital pairs with random orbits, following the expected dis- 
tributions (composite of uniformly random ascending nodes 
and the expected distribution of inclinations, i.e. the arc- 
cosine of a uniform distribution), and calculate differential 
orientations as above. The randomly generated distribution 
is then compared to the observed distribution through a 
Kolmogorov-Smirnov test. The probability that the popula- 
tions are drawn from the same distribution is 69%. The pop- 
ulation must therefore be described as being fully consistent 
with being random. Although there could hypothetically be 
some weak preference towards coplanarity, it is clearly far 
too small to make inferences about unseen inclinations in 
the systems. 

Since the dynamical history of multiple systems might 
be more violent than for an ordinary binary, it cannot be 
categorically excluded that a binary orbital inclination could 
be a better probe for the rotational plane of each single 
stellar component. However, it is clear that the viability of 
such a method would have to be carefully proven before it 
could be applied in practice. 

3.2 Large-number methods 

The use of asteroseismology to determine inclinations for 
a large number of targets in the context of transit proba- 
bilities is discussed in Beatty & Seager (2010), particularly 
with regards to SONG (Grundahl et al. 2008), with an es- 
timate that inclination determinations would require a time 
investment of approximately 10 nights per star. SONG is 
a planned network of 1-m class telescopes spread over all 
longitudes, equipped with high-resolution spectrographs for 
asteroseismology studies. For the transit survey envisioned 
in Beatty & Seager (2010) with 830 stars, the total require- 



ment becomes ~23 years, which seems excessive. However, 
in the context of the imaging survey studied here involving 
164 stars, we note that the requirement comes down to 4.5 
years, which would be well within the timescale for a launch 
of a dedicated Earth-like planet search mission, and hence 
feasible. Therefore, if asteroseismology can indeed provide 
inclinations to a reasonable precision, it appears a network 
like SONG could provide the required precursor science. 

A different possibility to determining inclinations for 
the bulk of the potential targets comes with the Space Inter- 
ferometer (SIM; e.g. Unwin et al. 2008). SIM would not only 
provide the inclinations and even phase angles, but also iden- 
tify which systems actually have planets in the first place, 
which is particularly useful if the frequency of Earth-like 
planets is low. On this note, it should be pointed out that 
SIM is in fact blind to the types of planets we have assumed 
for our simulations, since they would have a 1-year period, 
which is fundamentally indistinguishable from parallax er- 
rors. In any case, SIM would be hugely useful in determining 
inclinations for all stars that have any sufficiently massive 
planets. Taking the argument one step further, however, it 
might even be argued that since a main point of SIM would 
be to detect planets that are suitable for further study with 
direct imaging and spectroscopy, the target priorities of SIM 
should by themselves be partly influenced by a priori deter- 
minations of the inclinations, since those systems that are 
closest to face-on will be those that are easiest to study with 
imaging techniques. 



3.3 Additional utility of inclination 

In this paper, we have mainly focused on one strict util- 
ity from knowing planetary system inclinations a priori - 
maximizing the efficiency outcome from selection of targets 
for the detection of extrasolar Earth-like planets. However, 
knowing the inclination (and selecting for low-inclination 
systems) is useful for even wider purposes related to the 
study of such planets. With an a priori known inclination, it 
becomes easier to formulate a re-visit strategy once a given 
planet is detected, and for a close to face-on system, that 
schedule is minimally strict. Photometric or spectroscopic 
study of long-term atmospheric trends, in particular sea- 
sonal variations, benefits from low inclinations to maximize 
the continuous baseline. It might be argued that day-night 
variations would be harder to study in low-inclination sys- 
tems, if the obliquity is close to zero. However, Earth's obliq- 
uity is significantly non-zero, as are the Solar system planets 
in general, so this is unlikely to be a limiting factor. In ad- 
dition, interaction of an Earth-like planet with an exozodi- 
acal disk is likely to cause disk structures, which depend on 
the dynamical mass of the planet, thus potentially provid- 
ing a mass estimate independently of astrometry or radial 
velocity. Such structure is most easy to characterize at low 
inclinations, since ambiguity along the line of sight starts 
to occur at high inclinations. Finally, knowing the inclina- 
tion has positive spin-off effects on other areas of exoplanet 
studies such as for determining transit probabilities as in 
Beatty & Seager (2010), or for yielding real mass estimates 
for radial velocity planets. 
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4 CONCLUSIONS 

We have studied the impact of knowing the inclinations 
and/or exozodiacal brightnesses for nearby stars on the ef- 
ficiency of detection of extrasolar Earth-like planets. It was 
found that for a 60 mas IWA and a selection of 60 out of 
the 164 most nearby Sun-like stars, the efhciency can be 
boosted by 23% by knowing the inclinations of all the sys- 
tems, even if only >1000 zodi disks are known (as is the 
present status). If all >10 zodi disks are known but no incli- 
nations, the boost is 12%. A combination of the two provides 
a 32% improvement. These improvements may be smaller 
or larger depending on the exact distribution of exozodiacal 
brightnesses, and some mission parameters. Simultaneously 
determining all >10 zodi disks and all inclinations gives a 
near-ideal target selection. Attaining such a state of knowl- 
edge requires substantial effort but is probably realistic, and 
should be seriously considered as precursor science for ded- 
icated direct imaging missions. 

Although V sin i measurements and circumstellar mate- 
rial can provide inclination constraints for a fraction of the 
systems, it is probably necessary to use asteroseismology or 
space-based astrometry to provide inclination measurements 
for the bulk of nearby Sun-like stars. The orbital planes in 
multiple stellar systems is however not an acceptable proxy 
for the rotational planes of individual stars or the orbital 
planes of their planets. 
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